Light nuclei without a core 
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We introduce the no-core full configuration (NCFC) approach and present results for *He, ^^C 
and with the realistic NN interaction, JISP16. We obtain ground state energies and their 
uncertainties through exponential extrapolations that we demonstrate are reliable in *He where 
fully converged results are obtained. We find ^^C is overbound by 1.7 MeV and we predict the 
yet-to-be-measured binding energy of to be 68.1 ± 2.7 MeV. The extrapolated spectrum of ^''P 
is in reasonable agreement with known features of the ^''B spectrum. 

PACS numbers: 21.60.Dc, 21.60.Cs, 21.45.-v, 21.30.-x, 27.20.-|-n, 27.10.+h, 21.10.-k, 21.10.Dr 



INTRODUCTION AND MOTIVATION 

The rapid development of ab-initio methods for solving 
finite nuclei has opened the range of nuclear phenomena 
that can be evaluated to high precision using realistic 
nucleon-nucleon (NN) and threc-nucleon (NNN) interac- 
tions, even interactions tied to QCD [1, 2] where renor- 
malization is necessary [3] . Here we present methods for 
the direct solution of the nuclear many-body problem 
by diagonalization in a sufficiently large basis space that 
converged binding energies are accessed — either directly 
or by simple extrapolation. We do not invoke renormal- 
ization. We choose a harmonic oscillator (HO) basis with 
two basis parameters, the HO energy Ml and the many- 
body basis space cutoff A^maxj defined below. We assess 
convergence in this 2D parameter space. 

Such a direct approach may be referred to as a "No- 
Core Full Configuration" (NCFC) method. Given the 
rapid advances in numerical algorithms and computers, 
as well as the development of realistic non-local NN inter- 
actions that facilitate convergence, we are able to achieve 
converged results, either directly or through extrapola- 
tion. That is, we do not need to soften the NN interac- 
tion by treating it with an effective interaction formalism. 
Renormalization formalisms necessarily generate many- 
body interactions that significantly complicate the theory 
and are often truncated for that reason. Renormalization 
without retaining the effective many-body potentials also 
abandons the variational upper bound characteristic that 
we prefer to retain. While we omit all NNN interactions 
and higher-body forces for the present time wc do antic- 
ipate their need for a detailed understanding of the full 
range of experimental data. 

We adopt an m-scheme basis approach where the 
many-body basis states are limited by the imposed sym- 
metries — parity and total angular momentum projec- 
tion M, as well as by the cutoff in the total oscillator 
quanta above the minimum for that nucleus (iVmax)- In 
natural parity cases, M = (or i) enables the simulta- 
neous calculation of the entire spectrum for that parity 
and A^max- In many light nuclei, we obtain results for 



the first few increments of A^max and extrapolate calcu- 
lated observables from a sequence of results obtained with 
these A^max values. 

In our NCFC approach, the input Hamiltonian is inde- 
pendent of A^max; the computer requirements for a Con- 
figuration Interaction (CI) calculation are the same as 
that of the ab-initio No Core Shell Model (NCSM) [4] 
with a 2-body Hamiltonian rcnormalizcd to the chosen 
A'max basis. The NCFC results are obtained by taking 
the limit of A^max ^ oo. Both the NCFC and NCSM 
approaches guarantee that all observables are obtained 
free of contamination from spurious center-of-mass mo- 
tion effects. 



HAMILTONIAN INGREDIENTS, BASIS 
SELECTION AND METHOD OF SOLUTION 



In order to carry out the NCFC calculations, we 
require a realistic NN interaction that is sufficiently 
weak at high momentum transfers that we can obtain 
a reasonable convergence trend. The conventional Lee- 
Suzuki-Okamoto renormalization procedure of the ab- 
initio NCSM [4] develops effective interactions that pro- 
vide answers close to experimental observations. How- 
ever, the convergence trend of the effective interaction 
sequences with increasing A^max is not uniform and leads 
to challenges for extrapolation. Therefore, we select the 
realistic NN interaction, JISP16, that produces spectra 
and other observables in light nuclei that are in reason- 
able accord with experiment [5] . We include the Coulomb 
interaction between the protons. 

We cast the CI many-body problem with the "bare" 
interaction in the same harmonic oscillator (HO) basis 
and with the same definition of the cutoff as the ab-initio 
NCSM [4]. That is, the CI finite matrix problem is de- 
fined by A^max, the maximum number of oscillator quanta 
shared by all nuclcons above the lowest HO configuration 
for the chosen nucleus. The exact NCFC answer emerges 
in the limit A^max oo. This definition of the cutoff al- 
lows us to retain the same treatment of the ccnter-of-mass 
(CM) constraint that eliminates spurious CM excitations 
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as in the ab-initio NCSM. The CI matrix also depends 
on the HO energy, hfl. 

The CI approach satisfies the variational principle and 
guarantees uniform convergence from above the exact 
eigenenergy with increasing iVmax- That is, the CI re- 
sults for the energy of lowest state of each spin and par- 
ity, at any A^max truncation, are upper bounds on the 
exact NCFC converged answers and the convergence is 
monotonic with increasing iVmax- Our goal is to achieve 
independence of the two parameters as that is a signal for 
convergence — the result that would be obtained from 
solving the same problem in a complete basis. 

We employ the code "Many Fermion Dynamics — nu- 
clear" (MFDn) [6] that evaluates the many-body Hamil- 
tonian and obtains the low-lying eigenvalues and eigen- 
vectors using the Lanczos algorithm. 

By investigating the calculated binding energies of 
many light nuclei as a function of the two basis space 
parameters, we determined that, once we exclude the 
-^max = result, the calculated points represent an expo- 
nential convergence pattern. Therefore, we fit an expo- 
nential plus constant to each set of results as a function 
of A^max, excluding N^^nx = at fixed HQ, using the re- 
lation: 

£^gs(iVraax) = aexp(-c7V,„ax) + -Egs(oo). (1) 

EXTRAPOLATING THE GROUND STATE 
ENERGY — NCFC TEST CASES WITH *He 

We now investigate the convergence rate for the ground 
state energy as a function of iV^ax and hfl for "^He where 
we also achieve nearly exact results by direct diagonaliza- 
tion for comparison. In particular, we present the results 
and extrapolation analyses for *He in Figs. 1 through 3. 

The sequence of curves in Fig. 1 for ^He illustrates the 
trends we encounter in CI calculations when evaluating 
the ground state energy with the "bare" JISP16 interac- 
tion. Our purpose with ^He is only to illustrate conver- 
gence trends. The iVmax = 18 curve reaches to within 3 
keV of the exact answer that agrees with experiment. 

Next, we use these *He results to test our "extrapo- 
lation method A" as illustrated in Fig. 2. For extrapo- 
lation A, we will fit only four calculated points at each 
value of hil. However, in Fig. 2 we demonstrate the ex- 
ponential behavior over the range A^max = 2—16. Later, 
we will introduce a variant, "extrapolation method B" 
in which we use only three successive points for the fit. 
For extrapolation A, we select the values of to in- 
clude in the analysis by first taking the value at which 
the minimum (with respect to H^) occurs along the high- 
est A^max curve included in the fit, then taking one Ml 
value lower by 5 MeV and three hfl values higher by suc- 
cessive increments of 5 MeV. For heavier systems we take 
this increment to be 2.5 MeV. Since the minimum occurs 
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FIG. 1: (Color online) Calculated ground state energy of *He 
as function of the oscillator energy, hQ, for a sequence of A''niax 
values. The curve closest to experiment corresponds to the 
value A''niax ~ 16 and successively higher curves are obtained 
with A'max decreased by 2 units for each curve. 
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FIG. 2: (Color online) Calculated ground state energy of *He 
for A^max = 2—16 for JISP16 at selected values of hi}. Each 
sot of points at fixed hil is fitted by Eq. (1) producing the 
solid curves. Note the expanded energy scale. Each point 
is a true upper bound to the exact answer. The asymptotes 
Egs{oo) are the same to within 35 keV of their average value 
and they span the experimental ground state energy. 

along the A^^ax = 16 curve at hil = 20 MeV as shown in 
Fig. 1 , this produces the 5 curves spanning a range of 20 
MeV in hfl shown in Fig. 2. 

We recognize that this window of results in hQ values is 
arbitrary. Our only assurance is that it seems to provide 
a consistent set of extrapolations in the nuclei examined 
up to the present time. 

For the resulting 5 cases shown in Fig. 2, we employ an 
independent exponential plus constant for each sequence. 
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FIG. 3: (Color online) Extracted asymptotes and upper 
bounds as functions of the largest value of A^max in each set 
of points used in the extrapolation. Four (three) successive 
points in A'max are used for the extrapolation A (B). Uncer- 
tainties are determined as described in the text. Note the 
expanded scale and the consistency of the asymptotes as they 
fall well within their uncertainty ranges along the path of a 
converging sequence. 



perform a linear regression for each sequence at fixed Ml, 
and observe a small spread in the extrapolants that is 
indicative of the uncertainty in this method. Note that 
the results in Fig. 2 arc obtained with equal weights for 
each of the points. 

For extrapolation A, we will fit sets of 4 successive 
points due to a desire to minimize the fluctuations due 
to certain "odd-even" effects. These effects may be inter- 
preted as sensitivity to incrementing the basis space with 
a single HO state at a time while including two successive 
basis states affords tradeoffs that yield a better balance 
in the phasing with the exact solution. 

Next, we consider what weight to assign to each cal- 
culated point. We argue that, as iV,nax increases, we are 
approaching the exact result from above with increasing 
precision. Hence, the importance of results grows with in- 
creasing TVmax and this should be reflected in the weights 
assigned to the calculated points used in the fitting proce- 
dure. With this in mind, we adopt the following strategy: 
define a chisquare function to be minimized and assign a 
"sigma" to each calculated result at iVmax that is based 
on the change in the calculated energy from the previous 
-^max value. To complete these sigma assignments, the 
sigma for the first point on the A^max curve is assigned 
a value three times the sigma calculated for the second 
point on the same Gxcd-hQ trajectory. 

As a final clement to our extrapolation A strategy, we 
invoke the minimization principle to argue that all curves 
of results at fixed hfl will approach the same exact an- 
swer from above. Thus all curves will have a common 
asymptote and we use that condition as a constraint on 



the chisquare minimization. 

When we use exponential fits constrained to have a 
common asymptote and uncertainties based on the local 
slope, we obtain curves close to those in Fig. 2. The 
differences are difficult to perceive in a graph so we omit 
presenting a separate figure for them in this case. It is 
noteworthy that the equal weighting of the linear regres- 
sion leads to a spread in the extrapolants that is modest. 

The sequence of asymptotes for the ^He ground state 
energy, obtained with extrapolation A, by using succes- 
sive sets of 4 points in iVmax and performing our con- 
strained fits to each such set of 4 points, is shown in Fig. 
3. We employ the independent fits similar to those in Fig. 
2 to define the uncertainty in our asymptotes. In particu- 
lar, we define our uncertainty, or estimate of the standard 
deviation for the constrained asymptote, as one-half the 
total spread in the asymptotes arising from the indepen- 
dent fits with equal weights for each of the 4 points. In 
some other nuclei, on rare occasions, we obtain an outlier 
when the linear regression produces a residual less than 
0.999 that we discard from the determination of the to- 
tal spread. Also, on rare occasions, the calculated upper 
uncertainty reaches above the calculated upper bound. 
When this happens, we reduce the upper uncertainty to 
the upper bound as it is a strict limit. 

One may worry that the resulting extrapolation tool 
contains several arbitrary aspects and we agree with that 
concern. Our only recourse is to cross-check these choices 
with solvable NCFC cases as we have done [7]. We seek 
consistency of the constrained extrapolations as gauged 
by the uncertainties estimated from the unconstrained 
extrapolations described above. Indeed, our results such 
as those shown in Fig. 3, demonstrate that consistency. 
The deviation of any specific constrained extrapolant 
from the result at the highest upper limit A^max appears 
well characterized by the assigned uncertainty. We have 
carried out, and will present elsewhere, a far more exten- 
sive set of tests of our extrapolation methods [7]. 

As we proceed to applications in heavier nuclei, we face 
the technical limitations of rapidly increasing basis space 
dimension. In some cases, only three points on the iVmax 
curves may be available so we introduce extrapolation 
B. Our extrapolation B procedure uses three successive 
points in A^max to determine the exponential plus con- 
stant. We search for the value of H^l where the extrapo- 
lation is most stable and assign the uncertainty to be the 
difference in the ground state energy of the highest two 
points in A^max- As expected, since extrapolation B uses 
less "data" to determine the asymptote, it will have the 
larger uncertainty. Again, we trim the upper uncertainty, 
when needed, to conform to the upper bound. 

We present the behavior of the asymptotes determined 
by extrapolations A and B in Fig. 3 along with the ex- 
perimental and upper bound energies. In this case the 
results are very rapidly convergent at many values of 
producing uncertainties that drop precipitously with in- 



TABLE I: Binding energies of several light nuclei from exper- 
iment and theory. The theoretical results are obtained with 
the JISP16 interaction in NCFC calculations as described in 
the text. The uncertainties in the rightmost digits of an ex- 
trapolation is quoted in parenthesis where available. 



Nucleus /property 


Exp 


Extrap (A) 


Extrap (B) 


*He |-E;(0+,0)| [MeV] 


28.296 


28.299(1) 


28.299(1) 


''He |S(0+,1) [MeV] 


29.269 


28.68(12) 


28.69(5) 


''Li \E{l+,0)\ [MeV] 


31.995 


31.43(12) 


31.45(5) 


'^He \E{0+,2)\ [MeV] 


31.408 


29.74(34) 


30.05(60) 


'"C \E{0+,0)\ [MeV] 


92.162 


93.9(1.1) 


95.1(2.7) 


'''O |£(0+,0)| [MeV] 


127.619 


143.5(1.0) 


150 (14) 


"■^F \E{2-,2)\ [MeV] 


? 


68.1(2.7) 


69.4 




Exp (-92.162) 



Extrap (-93.9) 
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FIG. 4; (Color online) Calculated ground state energy of 
as function of the oscillator energy, hfl, for selected values of 
A'^max. The curve closest to experiment corresponds to the 
value Amax = 8 and successively higher curves are obtained 
with Amax decreased by 2 units for each curve. 



creasing A^max as seen in the figure. We note that the 
uncertainties conservatively represent the spread in the 
asymptotes since all the extracted asymptotes are consis- 
tent with each other within the respective uncertainties. 
The largest A^max points define the results quoted in Ta- 
ble I, a ground state overbound by 3 ± 1 keV. 



EXTRAPOLATING THE GROUND STATE 
ENERGY — NCFC FOR ^^C AND ^*F 
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FIG. 5: (Color online) Calculated ground state energy of ^'^C 
for Amax ~ 2—8 at selected values of hQ as described in the 
text. For each HQ the data are fit to an exponential plus a 
constant, the asymptote. The figure displays the experimen- 
tal ground state energy and the common asymptote obtained 
in extrapolation A. 




hQ. (MeV) 

FIG. 6: (Color online) Calculated ground state energy of ^''F 
for Amax = 0—6 with "bare" (solid lines) and effective (dashed 
lines) JISP16 interaction as function of the oscillator energy 
hfl. Shaded area shows a confidence region of extrapolation 
A predictions, crosses depict predictions by extrapolation B. 



In our investigations of the lightest nuclei [7] we ob- 
serve a marked correlation between binding energy and 
convergence rate: the more deeply bound ground states 
exhibit greater independence of hft at fixed A^max- Our 
physical intuition supports this correlation since we know 
the asymptotic tails of the bound state wave functions 
fall more slowly as one approaches a threshold for disso- 
ciation. This same intuition tells us to expect Coulomb 
barriers and angular momenta to play significant roles in 
this correlation. 



We proceed to discuss the ^^C results by introducing 
Figs. 4 and 5. The ^^C nucleus is the first case for which 
we have only the extrapolation from the A^max = 2—8 
results since the A^max = 10 basis space, with a dimen- 
sion of 7,830,355,795, is beyond our present capabihties. 
Thus, in order to illustrate the details of our uncetain- 
ties, we depict in Fig. 5 the linear regression analyses 
of our results spanning the minimum in hfl obtained at 
^max = 8. Extrapolation A produces overbinding by 
about 1.7 MeV. 
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Our next example is ^^F, an exotic neutron-deficient 
nucleus, the first observation of which is expected in an 
experiment planned in the Cyclotron Institute at Texas 
A&M University. In this case, we only attain the results 
up through iVniax — 6 presented in Fig. 6. The A^max = 8 
basis includes about 2 billion states and we will perform 
the respective calculations in the near future. There- 
fore we cannot estimate the accuracy of extrapolation B 
shown in the figure by crosses for different Ml values. In 
the case of extrapolation A, we need to include in the 
analysis the A^max = results and we obtain the binding 
energy prediction of 68.1 ± 2.7 MeV (shaded area in Fig. 
6) which is seen to be in a good correspondence with 
extrapolation B. It is interesting that, contrary to our 
NCFC approach, the trend of the conventional effective 
interaction calculations of the binding energy is mislead- 
ing in this case: the minimum of the respective Ml de- 
pendence is seen from Fig. 6 to shift up with increasing 
A^max indicating the development of a shallow minimum 
at A'max = 6 around fi,f2 = 12.5 MeV; the ground state 
energy at this minimum is above the upper bound re- 
sulting from the variational principle and the A^max = 6 
calculations with the "bare" JISP16 interaction. 

We performed also calculations of the excited states 
in "F. The results obtained with hfl = 25 MeV in the 
range of A^max values of 0-6, are presented in Fig. 7. 
We performed also the extrapolation B for the energies 
of these states. The respective excitation energies, i.e. 
the differences between the extrapolated energies and the 
extrapolated ground state energy, are also shown in the 
figure. The ^*F spectrum is seen to be in a reasonable 
agreement with the spectrum of the mirror nucleus ^^B. 
However we should note here that the spin assignments 
of nearly all states in the ^''B spectrum are doubtful. 

CONCLUSIONS AND OUTLOOK 

We present in Table I a summary of the extrapolations 
performed with methods introduced here and compare 
them with the experimental results. In all cases, we 
used the calculated results to the maximum A^max 
available with the bare JISP16 interaction. Our overall 
conclusion is that these NCFC results demonstrate 
sufficient convergence achieved for ground state energies 
of light nuclei allowing extrapolations to the infinite 
basis limit and estimations of their uncertainties. These 
convergence properties are provided by the unique 
features of the JISP16 NN interaction. The convergence 
rate reflects the short range properties of the nuclear 
Hamiltonian. Fortunately, new renormalization schemes 
have been developed and applied that show promise for 
providing suitable nuclear Hamiltonians based on other 
interactions with good convergence properties within the 
NCFC method [8] . Additional work is needed to develop 
th e corresponding NNN interactions. Also, further work 



is in progress to extrapolate the RMS radii. 
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